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A theoretical study is presented of surface waves at a monomolecular surfactant film between an 
isotropic liquid and a nematic liquid crystal for the case when the surfactant film is in the isotropic 
two-dimensional fluid phase and induces homeotropic (normal to the interface) orientation of the 
nematic director. The dispersion relation for the surface waves is obtained, and different surface 
modes are analyzed with account being taken of the anchoring induced by the surfactant layer, the 
curvature energy of the interface, and the anisotropy of the viscoelastic coefficients. The dispersion 
laws for capillary and dilatational surface modes retain structure similar to that in isotropic systems, 
but involve anisotropic viscosity coefficients. Additional modes are related to relaxation of the 
nematic director field due to anchoring at the interface. The results can be used to determine 
different properties of nematic-surfactant-isotropic interfaces from experimental data on surface 
light scattering. 
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I. INTRODUCTION 

The presence of a surfactant film at a fluid-fluid inter- 
face alters the dynamics of the interface. This is mani- 
fested in behavior of the interfacial waves, induced either 
externally or by thermal fluctuations [ESQ- The in- 
terfacial dynamics can be probed by measuring the light 
scattered on such surface waves (see the review by Earn- 
shaw Q). The scattering of light on surface waves is 
a powerful tool for probing the properties of surfactant 
films at fluid interfaces [H, ID, |6| , and a variety of systems 
have been recently investigated using this method (e.g. 
refs @, H i, El ELE , see also the review by Cicuta 
and Hopkinson [lj]). 

Recently, the application of surfactant films to modify 
the interfacial properties has been extended to the sys- 
tems in which one of the fluids is in liquid-crystalline 
phase (e.g. liquid crystal colloids [H[). The presence 
of a liquid crystal as one of the fluids complicates the 
problem of probing the interfacial properties by studying 
the dynamics of the surface waves for the following rea- 
sons. Firstly, there are additional degrees of freedom in 
the bulk of the liquid crystal phase due to its anisotropy. 
Secondly, the interaction with the surfactant film is more 
complicated due to anisotropic anchoring. Finally, the 
surfactant film in the anisotropic field created by the 
neighboring liquid crystal can itself show anisotropic be- 
havior, even if it behaves as an two-dimensional isotropic 
fluid at the boundary between isotropic fluids. 

A promising new direction for chemical and biological 
sensing devices has recently emerged which utilizes the 
properties of surfactant films self-assembled on the inter- 
face between water and a nematic liquid crystal. The 
surfactant film induces preferred orientation of the ne- 
matic director [ll| [lj], [HI • The adsorption of chemical 
or biological molecules at such interface can then lead to 
reorientation of the nematic director, enabling detection 
by an imaging system [ijj H3, HH, [H, [23| . 

In these methods, easy detection is limited to the sys- 



tems in which adsorption changes anchoring properties of 
the interface with respect to the adjacent liquid crystal 
phase quite considerably. Namely, the equilibrium an- 
choring angle should change in magnitude. The range of 
application of these systems could be made significantly 
broader, however, if a method were used that was sensi- 
tive to changes in the anchoring properties of the inter- 
face that did not necessarily result in nematic director re- 
orientation. For ex ample, the anchoring orientation may 
remain unchanged [15l 1 1 Qt | , the adsorption only changing 
the strength of the anchoring. 

If a small amount of an analyte is present in the wa- 
ter it may be adsorbed at the surfactant layer, provided 
the surfactant molecules possess appropriate chemical 
properties. Generally, such adsorption will result in a 
change in the elastic and viscous properties of the inter- 
face. Hence sensitive experiments which are able to de- 
termine the interfacial properties will allow much more 
detailed experimental insight into the properties of the 
interaction between the surfactants and the analyte than 
has hitherto been available, and experimental study of 
surface waves is a possible technique for this purpose. 

The theoretical description of surface waves at inter- 
faces between nematic and isotropic liquids was made 
back in 1970s [H, H, HI]. The results demonstrated 
that the spectrum of surface waves has a more compli- 
cated structure than in the isotropic case, and allows the 
use surface scattering experiments to determine proper- 
ties of nematic interfaces [13, [H, [H, Hl| . Since then, 
several theoretical and experimental advances have been 
made, and presently these systems remain a subject of 
investigation [32l. l33l. l34l l35l [361] . 

The present paper presents a theoretical study of the 
dispersion of the surface waves at a monomolecular sur- 
factant film between an isotropic liquid (e.g. water) and 
a nematic liquid crystal. The main distinguishing features 
of such interfaces, are (i) the anchoring induced by the 
surfactant layer, (ii) the curvature energy of the inter- 
face, (iii) reduction of surface tension due to surfactant, 
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and (iv) the anisotropy of the surface viscoelastic coeffi- 
cients. We base our treatment on the mechanical model 
for anisotropic curved interfaces by Rey [33] , which takes 
into account anchoring and bending properties of the sur- 
factant. We consider the case of the insoluble surfactant 
film that is in its most symmetric phase (isotropic two- 
dimensional fluid), and induces homeotropic (normal to 
the surface) orientation of the director. 

The paper is organized as follows. The continuum 
model used in the rest of the paper is set up in Sec- 
tion [TTJ In Section IIIII the dispersion relation for surface 
waves is derived. In Section IIVI the numerical solution 
of the dispersion relation is solved with typical values 
of material parameters, and dispersion laws for different 
surface modes are analyzed in absence of the external 
magnetic field, and the influence of the magnetic field is 
discussed in Section |V] 

The explicit form of the dispersion relation is written 
in Appendix [D] 



II. THE MODEL 

In this section we formulate the model of the 
surfactant-laden interface between an isotropic liquid and 
a nematic liquid crystal, used in the present paper, and 
write down the governing equations. We base our treat- 
ment upon the models of the nematic-isotropic interface 
by Rey [13, 5l[ , and well known hydrodynamic descrip- 
tion of isotropic liquids [lij and nematic liquid crystals 

mm. 

We consider the case when the surfactant film induces 
homeotropic (normal to the surface) orientation of the 
nematic director, which is usually true in a range of the 
surfactant concentrations [HI, [lEj, This case is the 
simplest to analyze, and, at the same time, the most 
important for biosensing applications where the direct 
change in anchoring angle cannot be always observed. 

We include optional external magnetic field in our 
study and limit our analysis by considering the direction 
of the magnetic field that does not change equilibrium 
orientation of the nematic director. 

We assume that the system is far enough from any 
phase transitions both in the surfactant film [3^| and in 
the nematic phase [4(| • Thus we avoid complications re- 
lated to the fluctuations of the nematic and surfactant 
order parameters and the divergence of viscoelastic pa- 
rameters near phase transitions. 

The surfactant films can exhibit rich phase behavior 
[39l | . and the form of the surface stress tensor depends 
upon the symmetry of the interface. However, this does 
not normally influence much the dispersion laws of the 
surface modes compared to the isotropic case 0. In the 
present paper we assume that the surfactant film is in the 
most symmetric phase (isotropic two-dimensional fluid). 
Although the symmetry of the film should break in pres- 
ence of the adjacent liquid-crystalline bulk phase, the film 
remains isotropic in equilibrium if the anchoring of the 



nematic is homeotropic, and symmetry breaking can oc- 
cur only due to fluctuations of the director field. If we in- 
troduce the order parameter for the film, the correspond- 
ing anisotropic contributions to the interfacial stress ten- 
sor would be of higher order in the fluctuations of the 
dynamic variables than is required in our linearized treat- 
ment, so such contributions can be omitted. 

We consider a surfactant layer at an interface between 
nematic and isotropic liquids to be macroscopically in- 
finitely thin. We assume that the surfactant film is in- 
soluble and Newtonian. This means that the model is 
applicable to systems in which the interchange of surfac- 
tant molecules between the interface and adjacent bulk 
fluids is small, and the relaxation of the orientation of 
surfactant molecules is fast compared to relaxation of 
surface waves. We also assume heat diffusion to be suffi- 
ciently fast so that the system is in thermal equilibrium. 
We do not consider systems where other effects, such as 
polarity, are important. 

We shall choose coordinate system in such a way that 
the unperturbed interface lies at a plane z = 0, the half- 
space z < is occupied by the uniaxial nematic liquid 
crystal, and the half-space z > is filled by the isotropic 
liquid. Other details of the geometry used in the present 
paper are summarized in Appendix [A"l 

The central equations in the present Section are the 
conditions for the balance of forces (Eq. (JT|) and torques 
(Eq. (O at the interface. The explicit form of these equa- 
tions depends upon the chosen macroscopic model, and 
the rest of this Section is devoted to formulation of the 
model used in the present paper. 



A. Balance equations 

The interfacial force balance equation is the balance 
between the interfacial force and the bulk stress jump: 



Here 



0. 



(1) 



(2) 



is the force per unit area exerted by the interfacial stress 



(3) 



is the force per unit area exerted by the isotropic fluid, 



(4) 



is the force per unit area exerted by the nematic liq- 
uid crystal, the subscript s indicates that the bulk stress 
fields in the isotropic liquid, and in the nematic, S w , 
are evaluated at the interface, k is the unit vector normal 
to the interface and directed into the isotropic liquid. 
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The interfacial torque balance equation can be cast as 

(5) 



where T s is the interfacial torque arising due to surface 
interactions, T N is the torque exerted upon the interface 
by the adjacent nematic liquid crystal. 

The explicit model for surface and bulk stresses and 
torques that enter Eqs Jl} and (Eq. ([5]) is expanded in 
the remainder of this section. 



where the tangential surface molecular field is given by 

hf = -Is ~ = - Is • (V s • M) , (14) 

V s is surface gradient operator, 5/ 5k. denotes variational 
derivative with respect to k. The elastic contribution to 
surface torque was written as 



-iSe 



e : S Se + V • Cc 



(15) 



B. Surface elastic stress and torque 

In this and the following subsections we summarize 
the equations for the surface stress tensor £ and surface 
torque vector T s . We represent these quantities as a sum 
of corresponding non-dissipative (elastic) and dissipative 
(viscous) contributions: 



iSt: 



-iSu 



(6) 



(7) 



To describe the non-dissipative contributions in the sur- 
face stress tensor, X Se , and surface torque vector, T Se , 
we use the equilibrium model proposed by Rey [33], 
which is summarized below. 

Rey considered the interface with the Helmholtz free 
energy per unit mass !F S of the form 



T h =T b (/,k,b), 



(8) 



where p s is the surface mass density, b is the second 
fundamental tensor of the interface (see Appendix [A"|l . 
The corresponding differential was written as 



where 



C s = -M-e s 



(16) 



is the surface couple stress, e is the Levi-Civita tensor, 
and es = —Is x k is the surface alternator tensor. 



C. Surface viscous stress and torque 

The viscous properties of interfaces between an 
isotropic fluid and a nematic liquid crystal were consid- 
ered in detail by Rey [4lJ , and the results are summarized 
below. 

The forces and fluxes that contribute to the dissipation 
function R were identified as follows: 



R = Ef " : S s + Ef w : A s 
dri|| 
~dt 



hf" 



I, 



hf 



kk- 



dn± 
~~dt 



(17) 



where Sf v and Sf v are, correspondingly, symmetric and 
antisymmetric parts of the surface viscous stress tensor 
£ s ", hf" and hf" are the components of the surface vis- 
cous molecular field tangential and normal to the surface, 



dT s = ^^+il.rfk+— :dh 



(p S ) 



where 



,2 dT s 



k.b 



is the interfacial tension, 



/ St d F S \ 



9k ) p s, h 



(9) 



(10) 



(11) 



is the tangential component of the capillary vector (Is is 
the surface projector), and 



M 



P" 



dh J 



(12) 



is the bending moment tensor. The elastic surface stress 
tensor was found to be 



iSe 



al s M b + hf e k 



(13) 



V.v 6 • L 



(V s v*)- 



(18) 



is the surface rate-of-deformation tensor ((• • • ) T denotes 
the transposed tensor), 



A y = - 



V s v 5 • I s - I s • (V s v s ) " 



(19) 



is the surface vorticity tensor, v s is surface velocity, 



and 



dn± 
~dT 



v s ■ V,n,i 



V s n_L 



(20) 



(21) 



are the total time derivatives of the components = 
I s • n and rij^ = kk • n of the nematic director field n, 
tangential and normal to the surface, correspondingly. 

Generally, presence of the surfactant film at the inter- 
face complicates the form of the entropy production due 
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to additional internal degrees of freedom of the surfac- 
tant, and to the anisotropy of the adjacent nematic liq- 
uid. However, if the surfactant film that is in its isotropic 
liquid phase and favors homeotropic anchoring of the ne- 
matic, the resulting anisitropic terms in the entropy pro- 
duction introduce corrections to the hydrodynamic equa- 
tions of higher order than linear, and therefore can be 
neglected in the linearized treatment. Since this is the 
case we are considering, we shall adopt the form of the 
entropy production (fT7|) in our model and use the form 
of the viscous contribution to the surface stress tensor 
derived by Rey [4l| , which is given by 

£ s * = afS s :n||n||n||n|| 

+ afnuN 5 + afN s n|| + af S s 

+ ofnpll -S s + alS s -n [m 

+ afn||n, t (n| r N 5 )+/5fl s (l s :S s ) 

+ /3 2 s [n,|n||(l s :S s )+I s (n||n|| :S S )], (22) 

where N s is the surface Jaumann (corrotational) deriva- 
tive [44| of the tangential component of the director , 
and af_ 7 , /3f_ 2 are nine independent surface viscosity co- 
efficients. In the isotropic case n = 0, the expression for 
the surface viscous stress tensor reduces to the viscous 
stress tensor of Boussinesq-Schriven surface fluid P, |45j 
with the interfacial shear viscosity 77,, given by 

Vs = ^, (23) 
and dilatational viscosity zrj s given by 

C. = ^+J9f- (24) 

The surface viscous torque, corresponding to Eq. (fl"7[). 
is given by [4l[ 

T Sv = -n x h Sv , (25) 
where the surface viscous molecular field h Sv is 
h Sv = 7 f A • ri|| + 7l ^N s + afri|| (n||ri|| : A s ) 

+ :S *)+ 7 f ± kk.^. (26) 

The viscosity coefficients can be expressed in terms of 
quantities af . We shall need only the expression for the 
tangential rotational viscosity: 

7ii = of - "I- (27) 

D. Anchoring and curvature energies 

To calculate explicitly the interfacial tension a 
(Eq. (fl"Q|) ). the tangential component of the capillary vec- 
tor £11 (Eq. (TTTj) ), and the bending moment tensor M 



(Eq. (p~2|) ) . we need to know the dependence of the surface 
free energy J- s on the orientation of the interface given 
by unit normal vector k, and on its curvature described 
by second fundamental tensor b. For small deviations of 
k and b from equilibrium, we can expand the free energy 
in powers of these quantities and truncate the series. The 
result can be represented as 

F s (p s ,k,h) = F t s (p s )+F!(k)+F?(h), (28) 

each of the contribution described below. 

The contribution Tf corresponds to the surface ten- 
sion a of the equilibrium interface (flat interface, adjacent 
nematic director normal to the interface): 

P S ^t = o. (29) 

The anchoring contribution to the surface free energy 
density, J 7 **, describes the energetics of the preferred 
alignment direction of the nematic director relative to 
the interface. For the homeotropic equilibrium anchor- 
ing, it can be written in terms of ny as follows: 

P 3 ?! = \wn\ + o(n\). (30) 

Such expansion applied to the widely used Rapini- 
Papoular form of the anchoring free energy density [46j 

P S T RP = ^(n-kf, (31) 

shows that these definitions of the anchoring strength 
coefficient have opposite signs: 

W = -W RP . (32) 

We shall use W as the anchoring strength coefficient to 
ensure that it is positive in the case of the homeotropic 
anchoring being considered. 

The third contribution to the surface free energy den- 
sity, , is caused by finite interface thickness, and is 
related to the difference of the curvature of a surfactant 
film from the locally preferred (spontaneous) value. The 
widely used form of this contribution is the Helfrich cur- 
vature expansion [4^, EH 

p s Tf = ~2kH 2 + 2k(H - Hj 2 + R,K + o({Vk) 2 ) . (33) 

Here the geometry of the interface is described by the 
mean curvature H and the Gaussian curvature K, and 
the material parameters characterizing the interface are 
the bending rigidity k, the saddle-splay (or Gaussian) 
rigidity ft, and the spontaneous curvature H . The term 
—2kH 2 guarantees that the curvature energy of a flat 
interface (H — 1 K — 0) is zero. 

E. Surfactant concentration 

To complete the description of the interface, we need 
the continuity equation for the surfactant concentration 
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v. For insoluble surfactants, the continuity equation 
reads: 



— + vV s • v s = 0. 
at 



(34) 



We shall extend the description of the dependence of 
the interfacial tension upon the concentration of sur- 
factant, presented by Buzza 0], to other parameters 
characterizing the interface (surface tension cf, anchor- 
ing strength W, bending rigidity k, saddle-splay rigidity 
R, spontaneous curvature H, and surface viscosities af , 
(if). For small deviation 5v = v — vq of the surfactant 
concentration v from its equilibrium value vq, these co- 
efficients can be written in form 



r \ _ , . da 

a[u) = a{u ) + -^-Si/, 
ov 



(35) 



and similarly for other quantities. Casting surface ve- 
locity v s as the time derivative of the small surface dis- 
placement u, 



v s = 



du 

IE 1 



(36) 



we obtain from the continuity equation Eq. (|34[) that 



Sv 



-v V s ■ u. 



(37) 



This allows us to represent the material parameters of 
the interface as 



<r{v) = oo + £qV s • u, 



W{v) = Wo + WiV s ■ u, 



K(v) = Kq + KiVj • U, 



H{v) = H a + i?iV s • u. 



(38) 



(39) 



(40) 



(41) 



In these formulas Oq — <t(vq), Wq — W(y Q ), n = k{v ), 
Hq = H(isq) are, correspondingly, the interfacial tension, 
anchoring strength, bending rigidity, and spontaneous 
curvature in the unperturbed interface, sq = —v^dajdv 
is the static dilatational elasticity, W\ = —v^dW/dv, 
K\ = —vodn/di 1 , and H\ = —v^dH/du are coefficients 
in the first order term of the expansion of anchoring 
strength, bending rigidity, and spontaneous curvature in 
powers of (V s • u) . There are similar expansions for Gaus- 
sian rigidity R and surface viscosities af , (if . 



F. Magnetic field 

Magnetic field 7i in the isotropic and nematic regions 
satisfies Maxwell equations [25j, |49( 



V 7i = 0. 



(43) 



Neglecting magnetization of the interface, the boundary 
conditions read 



(44) 



k-(W + 4 7 rM)| JV = k-(?i + 47rAl)| I . (45) 

Here the magnetization of the isotropic liquid is 

M\ I = X I n\ I , (46) 

where x 1 is the magnetic permeability of the isotropic 
liquid, the magnetization of the uniaxial nematic liquid 
crystal is (4fl | 



M\ N = X ±n + Xa(H- n)n\ 



N 



(47) 



where Xa is the difference of the longitudinal and 
transversal magnetic permeabilities of the nematic: 



Xa = X\\ 



(48) 



G. Isotropic liquid 



We assume both the isotropic liquid and the nematic 
liquid crystal are incompressible, so that their densities 
p 1 and p N , are constant. 

The linearized equations for the incompressible 
isotropic liquid are well known [42j ■ They are the conti- 
nuity equation 



V • v = 
and Navier-Stokes equations 
j9v 



Of 



V ■ S J , 



where the hydrodynamic stress tensor is given by 
S 7 = - P I + 2r/S, 



(49) 



(50) 



(51) 



where r\ is the shear viscosity of the isotropic liquid, I is 
the unit tensor, 



s = i 

2 



Vv+ (Vv) J 



(52) 



is the strain rate tensor. 

We assume the non-slip boundary condition for the 
velocities of bulk fluids adjacent to the interface, which 
means the equality of the velocity of surfactant, v s , and 
that of the bulk fluids at an interface, v|„: 



v x n = o, 



(42) 



(53) 
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H. Nematic liquid crystal 



III. DISPERSION RELATION 



To describe the dynamics of the nematic liquid crys- 
tal that is far from the isotropic-nematic transition and 
has small deviations from its equilibrium state, we shall 
use the linearized form of the Eriksen-Leslie theory 
40, 43, 50]. The linearized equations for the incompress- 
ible nematic liquid crystal are the continuity equation 
(|49p. the equation for the velocity 



s 



N 



and the equation for the director 

dSn 
~dt 

Here 



-7— = n • A + A (I — n n ) ■ S ■ n 



7i 



2 



Vv- (Vv) J 



(54) 



(55) 



(56) 



is the antisymmetric vorticity tensor, A is the reactive 
material parameter, 71 is orientational viscosity, h is the 
molecular field which, assuming Frank form of the elastic 
free energy of a nematic liquid crystal in magnetic field 



— [n • (V x n)] 



^3 r Ni2 

— [n x (V x n)] 



1 



Xa(n-H) 



has the linearized form 

h = (I n n ) ■ h*, 

where 

h* = KiVV ■ Sn - K 2 V x {n [n ■ (V x Sn)}} 



(57) 



(58) 



(59) 



+K 3 X7 x {n x [n x (V x Sn)]} + Xa (n • H)H, (60) 

K±, K2, and K3 are the splay, twist, and bend Frank 
elastic constants, correspondingly. 

The stress tensor can be represented as a sum of reac- 
tive and viscous (dissipative) contributions, 



Nr 



The linearized form of the reactive part is 



-pi + i (n h - hn ) - ^ (n h 



(61) 



hn ) . (62) 



The linearized viscous stress tensor of incompressible ne- 
matic is 



X Nv = 2^ 2 S + 2 (1/3 - v 2 ) (n S ■ n + n • Sn ) 

+ 2 (yi + v 2 - 2^3) n n n no : S. (63) 

The quantities V\, 71, and A can be expressed 

through more commonly used Leslie viscosity coefficients 
[50I l5l| . Note that equating r\ = v-y = v^ = v 3 recovers 
the viscous stress tensor 277S of the isotropic incompress- 
ible fluid (last term in Eq. (|5Tjl ). 



The aim of this Section is to construct the dispersion 
relation for the surface waves on the basis of the model 
set up above. We consider a surface wave with frequency 
u) and wavevector q = (q, 0, 0) propagating along x axis, 
and solve force balance equation, Eq. |T]), and torque 
balance equation, Eq. {5} using linearized form of the 
hydrodynamic equations written in Section ILT1 

In order to linearize the hydrodynamic equations, we 
represent pressure p — p(r, t) and the nematic director 
n = n(r, t), where r = (x, y, z) is the position in space, t 
is time, in form 



P = Pa + Sp, 



n = n + Sn, 



(64) 



(65) 



where Sp and Sn are the deviations of pressure and direc- 
tor from their equilibrium values po and n , correspond- 
ingly. The velocity v = v(r, t) is itself the deviation from 
zero equilibrium velocity. Homeotropic anchoring corre- 
sponds to 



n = (0,0,1). 



(66) 



For small deviations from the equilibrium, we shall use 
the hydrodynamic equations linearized in v, Sp, and <5n. 
We shall assume these quantities to be independent of the 
coordinate y (d y = d/dy = 0) and vanish at z — > ±00. 

The magnetic field can be also represented as Ti. = 
Hq + SH, where Hq = (0, 0, Hq) is the equilibrium value, 
and the deviation S7i can be found from the linearized 
form of the Maxwell equations (|42]) . (|43)l . The terms in 
the final equations, containing S7i, are of higher order 
than linear, so we shall use only the equilibrium value, 
and skip the '0' subscript, so that Ti. — (0,0, 

Substituting the interfacial free energy density (|28|) 
into Eqs [flO]). (fTTjl . (fl2]l. and ([14]). we find the contri- 
butions up to the first order in u (and its derivatives) 
and rij into surface tension 

a = a + ynjj-2Ki? 2 + 2k(H - H) 2 + RK, (67) 

bending moment tensor 

M = 2 [(« + «)# - kH] I s - fib, (68) 

tangential component of the capillary vector 

^=-Wn h (69) 

and tangential surface molecular field 

hf v = Wn {l - 2V S [(k + k)H-kH]+V s - (fib) . (70) 

The non-vanishing components of the surface viscous 
stress tensor () are 



2 Vs S s + (Cs - Vs) (Is : S s ) I s , (71) 
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The total interfacial force F s can be found by substitut- 
ing Eqs ©, JT3J), and 1(57 ]) -([71 ]) into Eq. ©, and has 
components 



where 



F x s = e d 2 x u x + ( Vs +( s )dlv x , 

Fy = Vsd 2 x v y , 

= (a + Wo)d 2 x u z + W d x Sn x 

+ ipd x ii x - K d x u z , 



(72) 
(73) 

(74) 



u = v/iu), we obtain balance equations for the force com- 
ponents in form 

- e o9 2 «f - ivv ( mI + l) v x - ur) (m 1 - q) v z - 



3 r 



-;/ 3 



. 1 + A 

-IUJ — 



E 

2 



N\ 



N\\ 



K 3 (mp) -K iq 2 + Xa H 2 



c 



0, (85) 



ip = 2 ( Ko iJi + mHo) 



(75) 



To write the explicit form of the force balance equations 
(HJ), we also need the expressions for the components of 
the force ([3|) exerted by the isotropic fluid, 



F x = v(dxV z +d z v x ) z=+0 , 

Fy = V(dzVy) z=+Q , 

Fl = {2 V d z v z -p) z=+0 , 



(76) 
(77) 
(78) 



and the components of the force (J2J) exerted by the ne- 
matic liquid crystal, 



F 



N 



1 + A 
2 

1 + A 



h x - v z (d x v z + d z v x ) 
h 



z=-0 



2 -y - V-sd z Vy 



(p - 2i/i5 z w z ) 2= _ . 



z=-0 



(79) 

(80) 
(81) 



The hydrodynamic fields v, p, n in the bulk isotropic 
and nematic liquids are found by solution of the hydro- 
dynamic expressions. The explicit formulas are presented 
in Appendices [B] and [Cl 

Next we introduce Fourier transforms in the x coordi- 
nate and in time as 



v(r,t) 



1 



dq / due tuJt - tqx w(q,z,u), (82) 



— OO — GO 



oo oc 



Sp (r, t) 



dq / duje wt -* qx p(q,z,uj), (83) 



-OO — OO 



OO OO 



n(r,i) = -i- 5 f dq f dtoe^-^h^z^) (84) 

(27T) J J 



(for brevity we shall henceforth omit arguments of the 
transformed functions). Performing Fourier-transform 
of the force balance equation ([T]), and substituting 



r] s q 2 Vy - r)m I; Uy - U S 



1 + A 



i=i 



E [*> (r 

»=i 



.JVJ 



W + XaH 2 



= 0, (86) 



„4-S 



~s 



w p + m w z 

777/ — q 



3 



J c 2 ^(A i -2^f")a JV|1 

1=1 



0,(87) 



where = eo + iuj (rj s + Cs) is the complex dilatational 
modulus, m 1 is defined in Appendix [Bl bv Eq. jB9j, and 

the quantities mf 11 , mf^, Cf 11 , Cf x , sf 11 , Bf x , and 
A l are define d i n Ap pendix Qby Eqs (JC52J), (jUTT|) . (jC29)l . 
(IU281) . (IC2T1) . (ICTOl) . and (IC20l) . correspondingly. 

To write the interfacial torque balance equation 
we cast the torque exerted upon the interface by the 
nematic liquid crystal, T ff , and the interfacial torque 
arising due to surface interactions, T s , entering the in- 
terfacial torque balance equation (O, in form 



and 



T w = kxh 



T s = k x h s , 



where the molecular field from the bulk 



-k- 



0(Vn) 



has linearized components 



.N 



-K 3 d z Sn s x , 



= -K 3 djn^ 



(89) 



(90) 



(91) 



(92) 



and the surface molecular field h s can be represented as 
a sum of elastic (h Se ) and viscous (h Sv ) contributions 



(93) 
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given by Eqs (fl4j) and (|26|) , correspondingly, and can be 
represented in components as 



h s x e = Wo (Sn s x + d x u z ) - KQ d 3 x u z + 4>d* 



h s y e = W 8n s y , 



lSv s & {Sn x + d x u z ) 

K =7l|| g t : 

Sv _ S dSn 



at ' 



(94) 
(95) 

(96) 
(97) 



Then the surface torque balance equations can be written 
as 



K 3 d z 8nl + Wo (5n x + - K d*u z + 

+ V<9>x + s: = 0, (98) 

d8n% 



Of 



- K z d z bn S y + W 8n y + 7 fj, = 0, (99) 

or, substituting the expressions (|C17|) . (|C19|) and ()C9|) . 

£ [(*amf " - Wo - iurf ) fl| + 



+ - ( VK + - ^m; v 11 +icj 7 f ) q" = 0(100) 



E ( W J 



Wo - ioj-fi ) b^c? 1 - = 0. (101) 



The interfacial force balance equations (Eqs (|85[) - (|57|) ) 
and the interfacial torque balance equations (Eqs (|100[) 
(fTOTj) ) form, with account of Eqs |C^8|) and (jC2ljj) . a ho- 
mogeneous system of linear algebraic equations in v x , v y , 
£>f , n x and . The dispersion relation is obtained from 
the condition of existence of a solution to these equations, 
i.e. the requirement for the determinant D (w, q) of the 
matrix of coefficients for this system to be zero 

D(u,q) = 0. (102) 

The equations (f56"| and (|101[) in v y and h y decouple 
from the equations ((551), ([57]) and (|100|) in uf , vf and nf . 
Therefore, the matrix of coefficients is block-diagonal, 
and the dispersion relation (Eq. (|102p ) is equivalent to a 
pair of relations for (x, z) and y directions: 



£>H (w,?) = 0, 
D ± (w,g) = 0, 



(103) 
(104) 



where Z?|| (w,<?) is the determinant of the 3x3 matrix 
My of coefficients for the equations (|55|). 1(57) and ([TOO)) , 
and {^,q) is the determinant of the 2x2 matrix of 
coefficients for the equations (l86l) and (jlOip . 

The explicit form of the dispersion relations is pre- 
sented in Appendix [D] and can be readily used for the 
numerical analysis of surface modes. 



IV. SURFACE MODES 

In this Section the dispersion equation, which is pre- 
sented in Appendix [D] is solved numerically, and surface 
modes of different types are analyzed. For simplicity, 
we assume the density of the isotropic liquid, p 1 , to be 
small enough to be neglected (e.g. nematic-surfactant- 
air interface). We also assume that the magnetic field is 
absent. 

The surface modes can be easily classified at low 
wavevectors q. Expansion of the dispersion relation in 
powers of the wavevector q is a straightforward exercise 
in algebra, and the resulting modes are described below. 

Firstly, there is a transverse capillary mode, which 
has the dispersion law similar to that in the case of an 
isotropic liquid-liquid interface 0, H[ : 



uc (<?) 



'er <r 



P 



N 



o\q 



3/2 



(105) 



The principal contribution to this mode at large wave- 
lengths arises due to the restoring influence of surface 
tension do, and the predominant motion is in the direc- 
tion normal to the interface (z). The differences from the 
isotropic case, related to anisotropy of viscous dissipation 
in the nematic, appear in higher orders in q. 

The dilatational (or compressional) mode with pre- 
dominant motion in the direction along wave propaga- 
tion (x) arises in presence of surfactant layer due to the 
restoring force provided by the dilatational elastic mod- 
ulus eo. The dispersion law for this mode can be written 
as 



Ud (<?) 



■ ? A 



.,JV 



1 1/3 



W J 



■oiq 



, 4 /3 



(106) 



where the Miesowicz viscosity r/Q is given by [5 



% M =^3 



a+Ar 



-71. 



(107) 



The difference from the dispersion law for the dilatational 
mode in the case of a surfactant film at the interface 
between isotropic fluids, given by 0, [f| 

J2„4\ 1/3 



arises due to anisotropy of viscous dissipation in nematic. 

A new mode, specific to the nematic, is driven by relax- 
ation of the director field to equilibrium due to anchoring 
at the interface and has the disperion law 



wjv (<?) = — g- + o(q ) 
7i 



(109) 



Such relaxation is present even in absence of motion of 
the interface (e.g. when the interface is solid), so that 
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iojn does not vanish at q — > 0. For nematic-isotropic 
interfaces, the corresponding motion of the interface is 
induced by backfiow effects. 

Finally, behavior of the in-plane shear mode, with mo- 
tion in y-direction, is also governed by relaxation of the 
nematic director due to anchoring. The corresponding 
dispersion law 



t \ lW ° , ^ on 
(q) = + o [q ) 

7i 



(110) 



appears to be different from the isotropic case, where 
the damping of the in-plane shear mode in absence of 
anchoring is governed by the surface viscosity 77^ Q. 

Gravity g, so far neglected in our analysis, becomes 
important at wavevectors 




(111) 



and can be taken into account by adding the hydrostatic 



pressure term —g\p 



N 



to Eq. (|D28| . which corre- 



sponds to the additional contribution — g |p JV — p 1 1 u z to 
the vertical component of the force, Eq. (f?4"|) . The re- 
sulting dispersion law for transversal mode is given by 
expression 



u G (9) = Vgq + olq 



,1/2 



(112) 



which describes well-known gravity waves [42| . 

In the opposite case of large wavevectors, the curvature 
energy becomes important. Analysis of Eqs (|74f and (|94)) 
yields the characteristic values of q 



Q ~ q K 




and 



Wn 



q ~ oty 



(113) 



(114) 



below which one can neglect in the dispersion relation the 
terms containing bending rigidity k and its derivatives 
with respect to surfactant concentration, given by ip. 

Usually <fy > q K , and the range of q in which both 
gravity and curvature contributions become small, given 
by 



q g < q < q K 



(115) 



is rather wide. For typical values \p N — p 1 ~ 10kg/m 3 , 
g ~ 10m/s 2 , cr ~ W ~ 10- 2 J/m 2 , n ~ 10- 18 J, the 
equation (|115p reads 1cm -1 <C q <C 10 6 cm _1 , which 
includes the range of wavevectors typically probed by 
surface light scattering experiments. 

To obtain the dispersion laws for surface modes at 
larger values of the wavevector q, the dispersion equa- 
tion must be solved numerically. The numerical solu- 
tion presented below uses the following typical values of 
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Figure 1: Dispersion law uo (q) for different surface modes 
in absence of gravity, obtained by solution of the dispersion 
relation (|103|) with the values of the parameters given in the 
text. Numbers 1, 2, 3 denote transverse, dilatational, and ne- 
matic director relaxation modes, correspondingly. Prime and 
double prime denote real (solid line) and imaginary (dashed 
line) parts of u), correspondingly. 



the material parameters when it is not indicated other- 
wise. For the nematic liquid crystal we use the param- 
eters of 4-n-pentyl-4'-cyanobiphenyl (5CB) at 26°C 
the density p N — 1021.5kg/m 3 , the elastic constants 
K x = 5.95 • 1(T 12 N, K 2 = 3.77 • 1(T 12 N, K 3 = 7.86 • 
1CP 12 N, the Leslie viscosities ot\ = —6.6 • 10~ 3 kg/(m-s), 
a 2 = -77.0 • 10- 3 kg/(m-s), a 3 = -4.2 • 10" 3 kg/(m-s), 
a 4 = 63.4 • 10- 3 kg/(m-s), a 5 = 62.4 ■ 10~ 3 kg/(m-s), 
«6 = — 18.4- 10 -3 kg/(m-s). The viscosity coefficient used 
in the present paper can be calculated from the Leslie 
equations [EH and equal v\ = 50.4 • 10~ 3 kg/(m-s), 
u 2 = 31.7 • 10" 3 kg/(m-s), v 3 = 19.96 • 10~ 3 kg/(m-s), 
71 = 72.8 • 10" 3 kg/(m-s), A = 1.115. We use the value 
of the bending rigidity kq = 10~ 19 J which is typical for 
surfactant layers J53|. For other parameters we use the 
following typical values: 7-? ~ 5 • 10~ 9 kg/s, Hq = 0, 
e = l(T 3 N/m, 7} s = lfr 8 N/m, <r = 5 ■ lO^J/m 2 , 
W = 20 ■ 10~ 3 J/m 2 . 

The dispersion law ui (q) for different surface modes in 
absence of gravity, obtained by solution of the dispersion 
relation (|103[) with the values of the parameters given 
above, is presented in Figure [TJ At low q the dispersion 
of for modes 1, 2, 3, as denoted Figure [IJ is in good 
agreement with approximate formulas (|105p . (|106p . and 
(|109p . correspondingly. The noticeable discrepancy in 
behavior of capillary and dilatational modes appears at 
q ~ 10 cm -1 , and the damping of surface waves becomes 
large at larger q, which is qualitatively similar to the case 
of the interface between isotropic liquids. 

The results presented in Figure Q] suggest that in the 
typical range of q probed by surface light scattering ex- 
periments (100 cm -1 < q < 2000 cm^ 1 ), the approximate 
expressions (|105l) . (|106p do not describe well the disper- 
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Figure 2: Dispersion law w (g) for different surface modes in 
presence of gravity (/ = 9.8 m/s 2 , obtained by solution of the 
dispersion relation (|103[) with the values of the parameters 
given in the text. Numbers 1, 2, 3 denote transverse, dilata- 
tional, and nematic director relaxation modes, correspond- 
ingly. Prime and double prime denote real (solid line) and 
imaginary (dashed line) parts of u>, correspondingly. Vertical 
dotted line corresponds to the value of q g given by Eq. (jlllfl . 
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Figure 3: Dependence of the real (solid line) and imaginary 
(dashed line) parts of the frequency of the mode 1 (as defined 
on Figure [1} upon the bending rigidity Ko, calculated at q = 
1000 cm -1 in absence of gravity. Vertical line corresponds to 
the value of kq that satisfies Eq. (|113|l . 



sion curves, and accurate solution of the dispersion equa- 
tion should be used instead. 

Figure [5] presents the dispersion law uj (q) for different 
surface modes obtained by solution of the dispersion re- 
lation (|103[) in presence of gravity g — 9.8 m/s 2 . In agree- 
ment with the discussion above, the influence of gravity 
on the dispersion laws is small at q 3> q g , where q g is 
given by Eq. Qlll[) . 

If the bending rigidity k is large, its influence becomes 



Figure 4: Dependence of the real (solid line) and imaginary 
(dashed line) parts of the frequency u of the mode 3 (as de- 
fined on Figure [TJ, normalized by Wo/yf (see Eq. ((T09J), 
upon the anchoring strength Wo, calculated at q = 100 cm -1 
in absence of gravity. 
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Figure 5: Dependence of the real (solid line) and imaginary 
(dashed line) parts of the frequency ui of the in-plane shear 
mode, normalized by Wo/ji (see Eq. qilOp ), upon the an- 
choring strength Wo, calculated at q = 100cm -1 in absence 
of gravity. 



noticeable, as it is demonstrated in Figure [3] For k ~ 
kT, typical for surfactant films, the value of q K , given 
by eq. (|113[> . corresponds to wavelength close to atomic 
scales, and curvature energy can be neglected in typical 
surface light scattering experiments, in agreement with 
the discussion above. 

The dispersion law for the modes governed by relax- 
ation of the nematic director field in x and y directions 
due to anchoring of the nematic director at the interface, 
obtained by numerical solution of the dispersion equation 
with the values of the parameters given above, are well 
described by the equations (|109p and (| 1 1 0[) . However, 
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director tends to be oriented along the field rather than 
to be advected with the nematic liquid. The practical use 
of this effect is, however, limited, because at short wave- 
lengths extremely large magnetic field is required, and at 
long wavelengths gravity becomes dominating (Eq. (|112p . 

In principle, magnetic field can also influence surface 
waves through change in the properties of the interface 
(e.g. surface tension) due to the magnetization of the 
surfactant. Separate study is required to estimate the 
magnitude of this effect. 




1 



VI. CONCLUSION 



Figure 6: Dependence of the real (solid lines) and imaginary 
(dashed lines) parts of the frequencies of the modes 1 and 2 
(as defined on Figure [TJ upon the magnetic field, calculated 
at q = 10 cm -1 in absence of gravity. Vertical dotted line 
corresponds to XaT~L 2 = w 7i ( see Eq. (|116|l ). 

as the anchoring strength becomes smaller, other mech- 
anisms start to take over, as demonstrated in Figures [4] 
and [5] 

V. INFLUENCE OF MAGNETIC FIELD 

In this Section we discuss how the surface modes de- 
scribed in Section HVl are altered in presence of the exter- 
nal magnetic field directed normally to the surface (along 
z axis). 

The external magnetic field effectively acts on the ne- 
matic molecules as an additional molecular field (see 
Eq. (|59|) ). and the primary counteracting mechanism is 
provided by orientational shear relaxation. Thus we may 
expect the influence of the magnetic field become notice- 
able at 

Xa H 2 ~ w lv (116) 

The results of the numerical solution of the dispersion 
equation in presence of magnetic field, presented in Fig- 
ure [6j confirm that noticeable change in dispersion of 
capillary and dilatational modes arises only around the 
value of the field given by Eq. (|1 16[) . The change due to 
magnetic field in modes governed by anchoring is found 
to be negligibly small. 

At low q the dispersion of a capillary mode in strong 
magnetic field is different from the law (|105p and is given 
by 

/ (<TQ + Wq) g 3 " , / 3/2 \ 

wc v — — +0 \ q )■ (117) 

The frequency of this mode becomes sensitive to the an- 
choring properties of the interface, because the nematic 



We have obtained the dispersion relation for the sur- 
face waves at a surfactant-laden nematic isotropic in- 
terface for the case when the surfactant film induces 
homeotropic (normal to the surface) orientation of the 
director, and the surfactant film is in the isotropic two- 
dimensional fluid phase. We have analyzed the dispersion 
law of different surface modes analytically in long wave- 
length limit, and numerically in broader range of wave 
vectors, using typical values of the material parameters. 

At long wavelengths the dispersion of capillary, dilata- 
tional (or compression), in-plane shear, and director re- 
laxation modes is described by equations (| 1 05[) (or (|112p h 
(fTOfJf . (fTTOf . and (fT09|) . correspondingly. At smaller 
wavelength, the solution of the full dispersion relation 
should be used. Gravity influences the transversal mode 
at small wavevectors (Eq. (|lllj) ). and curvature energy of 
surfactant can be neglected if wavevector is not too large 
(Eq. (|113p V For all modes, the influence of the external 
magnetic field directed normally to the interface is small. 

The influence of the magnetic field should be more pro- 
nounced if the direction of the field does not coincide with 
equilibrium nematic director. In this case the dispersion 
law for surface modes may be expected to be quantita- 
tively different due to anisotropy of viscous dissipation 
in nematic, and different anchoring energy. The results 
of the present paper can be readily extended to the case 
of arbitrary direction of the external field and to other 
types of nematic anchoring. 

Other possible developments, which may increase the 
range of accessible systems and conditions, is the exten- 
sion of the results to wider range of the states of the 
surfactant film, and the study of the effects which may 
be caused by the phase transitions in the surfactant film 
and bulk liquid crystal. 

Dependence of the dispersion waves upon the param- 
eters of the interface suggests the surface light scatter- 
ing on a surfactant-laden nematic-isotropic interface as 
a potential method for determining of the properties of 
surfactant-laden nematic-isotropic interfaces, and as a 
possible candidate for a chemical or biological sensing 
technique. 
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Appendix A: DIFFERENTIAL GEOMETRY OF 
THE INTERFACE 



The geometrical description we use is similar to that of 
that presented in works Q and [37J • We choose the plane 
z = to coincide with the unperturbed interface, the 
half-space z < to be occupied by the uniaxial nematic 
liquid crystal, and the half-space z > to be filled by the 
isotropic liquid. 

Let the position of a fluid particle at the interface be 
r = ro + u, where ro = (xo,yo:0) is its position on the 
undeformed interface [z = 0), and u = u(ro) is the dis- 
placement vector with components (u x ,u y ,u z ). We shall 
use Xq and yo as surface coordinates and denote them as 
s a , a and other Greek indices taking values 1 and 2. 

The position r of fluid particles at the interface in 3D 
space can be cast as 



(Al) 



The surface tangent base vectors a a = dr/ds a , corre- 
sponding to the chosen surface coordinates, can be writ- 
ten in terms of the components of the displacement vec- 
tors: 



Or 



ai = -r— j- = (1 + d x u x , d x u y , d x u z ) 



(A2) 



and 



a 2 = = (d y u x , 1 + dyUy, d y u z ) . (A3) 



ds 2 

The surface metric tensor 



a a j3 = a a • a/3 



has determinant 



1 + 2d x u x d x u y + d y u.j 

X Uy + dyU X 1 + 2dyUy 



0(u 2 



(A4) 



a = dct(a Q/3 ) = 1 + 2 (d x u x + d y u y ) + 0(u 2 ). (A5) 

The corresponding reciprocal base vectors a Q and metric 
tensor a a/3 take form 

ds 1 

a 1 = — = (1 - d x u x , -d y u x , d x u z ) + 0(u 2 ), (A6) 



ds 2 

a 2 = — = (-d x u y , 1 - dyUy, d y u z ) + O (u 2 ) , (A7) 



a af3 = a" •a' 3 



(A8) 



~0 X Uy dy1l x 1 Qidyliy 



+ 0(u 2 ). 



1 
1/' 



(A9) 



We write the unit vector k, normal to the interface and 
directed into the isotropic liquid, as 

k = ai X&2 . = (-d x u z , -d v u z , 1) + 0(u 2 ). (A10) 
|aixa 2 | v ' 

We shall also define the dyadic surface idem factor 



1 d x u z 



I s =a Q a Q =l 1 d y u z \+0(u 2 ). (All) 

\d x u z d y u z 

the surface gradient operator 
V s =I S -V 



d x + {d x u z ) d z 

dy + (dyU z )d Z ) +0(U 2 ). (A12) 

(d x u z ) d x + (d y u z ) d y . 



and the second fundamental tensor 

d 2 x u z d x d y u z 0\ 
b = -V s k = ( d x d y u z d 2 y u z + 0{u 2 ). (A13) 
0/ 

The mean curvature H and Gaussian curvature K are 
given by 

H = 1(1, : b) = 1 {d 2 x u z + d 2 y u z ) + 0(u 2 ), (A14) 



K = ^P e ^ baib0S = Q{u 2y (A15) 

Other useful identities include the surface projection 
nil of a nematic director field n, Eqs (f65|) and ([66]) . 

' '5n x + d x u z s 

nil = Is ■ n = I 5n y + d y u z I + o (Vu, 8n) , (A16) 



and its surface divergence 

V s ■ n|| = d x Sn x + d y 6n y + d 2 u z + d 2 u z + o (Vu, Sn) 



(A17) 



Appendix B: BULK SOLUTION FOR ISOTROPIC 
LIQUID 

This appendix presents the solution to the linearized 
hydrodynamic equations in bulk isotropic liquid, ob- 
tained by Kramer 

Substitution of Eq. ([82]) into Eq. gH) yields 



iqv x + d z v z = 0. 



(Bl) 
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Substituting Eqs (JUJ) and (83) into Eqs ((50J) (J52l) , we 
obtain 

[iujp 1 + V (q 2 -d 2 z )] v x = iqp, (B2) 

[tup 1 + V {q 2 ~ d 2 )} v y = 0, (B3) 

[iuj p 1 + r 1 (q 2 ~d 2 )]v z = -d z p, (B4) 



where equation (|B3[1 is decoupled from other equations. 
The general solution to Eqs (|B1|) - (|B4[) vanishing at z — * 
oo can be written as 



h y = -K 2 q 2 n y + K 3 d 2 z fi y + XaH 2 n y . (C2) 

Substituting them into Eqs (154")) , (julj) ([(IB")) , we obtain the 
following linear differential equations, 

[iujp N + (2v 2 - iAj) q 2 - v z d 2 z \ v x - 1 - - dji x = iqp, 

(C3) 



1 - A 



[iup N + v 2 q 2 - v z d 2 z \ v y - ^—^d z h y = 0, (C4) 



with 



C IJ -e- m z , 



syl\\ —qz , /^yi" —m 5 

v z = C x e w + U 2 e 



m 



(B5) 



(B6) 



(B7) 



, 2 + ^ 1/2 , Rem J >0. (B9) 
V J 



The quantities C-f", Cj", and C 7 ^ are functions of q and 
w and are determined by the boundary conditions at the 
interface as follows: 



m 1 — q 



(BIO) 



[iiop N + v z q 2 - (2vi - is 3 ) d 2 ] v z - iq l+ ^ h x = -d z p, 

(C5) 

which are analogous to Eqs (|B2|) (|B4|1 for isotropic liq- 
uids. Equation (|55|) for the director after Fourier trans- 
form gives two equations, 

1 + A„ _ 1 - A . _ l r 

iun x = —r—a z v x H —iqv z H /i^ (C6) 

2 2 71 



(B8) and 



1 + A 



iojn v 



d z v v 



1 r 



2 7i 



(C7) 



where /i^ and h y are given by Eqs (|C1[) and (|C2|) . 

Thus we have six linear differential equations (Eqs (|B1|) 
and (|C3p - (|C7p ) for six dynamic variables (pressure, three 
components of velocity, and two components of director) . 
Equations (|C4[) and (|C7[) for v y and n y decouple from the 
others, their general solution vanishing at z — ► — oo can 
be cast as 



iVX mf 1 ! 



(C8) 



(Bll) 



?ti j — q 



(B12) 



where the superscript S indicates that the values of the 
corresponding dynamic variables are taken at z — > 0. 



Appendix C: BULK SOLUTION FOR NEMATIC 
LIQUID CRYSTAL 

In this appendix the solution is presented to the lin- 
earized hydrodynamic equations in bulk nematic liquid 
crystal. 

For the equilibrium director along z axis (Eq. (|66p ) 
the Fourier-transform, similar to equations (|B"!2"l) - (f81)) . of 
the linearized molecular field (Eq. (|58[) ). h, has non-zero 
components 



K iq 2 fi x + K 3 d 2 z h x + xan 2 fi x , (ci) 



N.L m " z 

t" 1 



where 



(1 + A) 7im 



2 lW7l +% 2 ~if 3 (mf 1 ) - X aU 7 



(C9) 



(CIO) 



mf x =( M ^) 1/2 , Remf ± >0, (Cll) 



and /x^ 1 , i = 1,2, are the roots of the quadratic equation 



4 (I*- 1 ) + «fV + 4=0, 



where 



1-A 2 



-7i - f 3 # 3 , 



(C12) 



(C13) 



14 



aj 1 = (iuj! + K 2 q 2 - XaU 2 ) v 3 + {iup N + v 2 q 2 ) K 3 
1-A 2 



-7i (K 2 q 2 ~ XaU 2 ) , 



(C14) 



(icop N + v 2 q 2 ) (turfr + K 2 q 2 - Xa H 2 ) . (C15) 



The general solution to the equations (|B 1|) . (|C3|) . (IC5I) . 

and (1C6|) vanishing at z — -> — oo can be cast as 



of = -zw7i [iwp^ + 2 (^i + ^ 2 - ^3) o 
'1 + A 2 



-71 -2{V\+V2- 1/3) 



(if l0 2 - Xa^ 2 ) 



1-A 2 



-7i ~vz\q 



K 3 q 2 , (C26) 



e. = £c?V 



p=£ac, 



(C16) 



(C17) 



(C18) 



(C19) 



zjj = iLUjiq 2 (iuop N + v 3 q 2 ) + [iwp 
'1-A 2 



N 



-7i - ^3 9 



The quantities Cf 11 and Cf ^ are functions of q and 
a; and are determined by the boundary conditions at the 
interface as 



(K lQ 2 - Xa H 2 ) q\£27) 



Rt-B£ ' 



(C28) 



where 



At 



1 f 1 + A 



N\\ 



K lQ 2 - K 3 ( mf !t ) — XaH 2 



B 



iujp N + v 3 q 2 - (2i/i - t/ 3 ) (mf 



II _ *7i 



(1-A)g 2 -(1 + A) (mf 11 )' 



2Q 



^ 7l +^ l0 2 -X 3 (™f") -XaH 2 



1 /2 

(/if) , Remf'x), (C22) 



(C20) 



(C21) 



and fx 1 }, i = 1, 2, 3, are the roots of the cubic equation 
4 (jA +c& ^ 2 +44 + 4 =0, (C23) 



where 



1-A 2 



-7i - v 3 K, 



(C24) 



a 2 = iu)~{iv 3 



1-A 2 



71-^3 (K iq 2 - Xa n 2 



+ S UjJP 



N 



1 + A 2 



7i-2 (vi +v 2 - v 3 ) 



q 2 K^C2h) 



Cf" = [i q (Bl-Bl)vl+( 



N\\ N\\\ ~S 

m Q — rrir, nz. 



m N jBl-< ] Bl «f /A, 



(C29) 



where 



A = {B\-Bl)m^ + {B\ 



(B\-B^m»\ 



(C30) 



and expressions for C 2 ^, Cf ", and cf " are obtained 
from Eqs (|C28[) and (|C29|) by cyclic permutation of sub- 
script indices. 



Appendix D: EXPLICIT FORM OF DISPERSION 
RELATION 



To write the explicit form of the dispersion relations 
Pllg)) and Q3H> , we recast equations ([CM)) and ([C29]) in 
form 



Cf 11 = + L™v t , +L\ n x) n x (D2) 



(Dl) 
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where 



L 



r (vy 

r (vy 
J 2 

-(ny 

(ny 
2 

- (vx 
J l 

- (vx 
J 2 

(vx 
3 

(vz 



- (nx 



- (nx 



- (nx 



B 2 - 
Bj 

1 

Bf^Bj' 
1 

iq f Bj - B\ 



(Bl-Bl) 



"I 



h(b\-b^ 




A 


N\ 

m 3 


Bl-m^Bl 




A 


N\ 

m 1 


Bi-mfBf 






N\ 

m 2 


Bi-mfBl 




A 


JV| 


N\\ 

— m 3 




A 


JV| 

m 3 


AT|| 

— m 1 




A 


JV| 

m 1 


AT|| 

- m 2 



(D3) 
(D4) 
(D5) 
(D6) 

(D7) 

(D8) 

(D9) 
(D10) 
(Dll) 
(D12) 
(D13) 
(D14) 
(D15) 



A is given by Eq. (|C30|) . Bj- and Bj are given by 

Eqs (|C10|) and (|C21j) . mf ± and mf" are given by 
Eqs (|C22[) and (jClip . correspondingly. 

Then the dispersion relation (| 104[) can be written as 



det M 1 - = 0, 



(D16) 



where M 1 - is 2 x 2 matrix of coefficients for equations 
flgBJ) and (HUB) 



M^ M^ 



with the following components: 



M 



ii = -»M - ?7 m 



J +E 



,AT_L 



1 + A 



X 3 (mf ^) 2 - K 2 q* + Xa U 2 ] Bi) L^,(D18) 



= ^ (# 3 mf ^ - Wo - z^ 7 f) ^^r' , (D20) 



M 2 i = £ {K 3 m^ -Wo- iun?) B^L ( r\ (D21) 
»=i 



The dispersion relation (|103|) can be written as 

det Ml 1 =0, (D22) 

where M" is 3 x 3 matrix of coefficients for equations 
dHSD, 457), and (JTJHJ) 



Af 11 



m| 2 m| 3 
m^ m| 2 m| 3 
Va4 m| 2 m| 3 , 



(D23) 



with the following components: 

-Cq5 2 — iujr] (m 1 + q) 



Ml 



3 r 



-aE K") +^ 
»=i L 



r («a0 . . 1 + A 
L) ' + — - — x 



E 



^ 3 (mf") -K iq 2 + XaH 2 



Bl4 vx \ (D24) 



M-[' 2 = —u)i] {m 1 — qj 



L) + iuj- 



E K") + 



1 + A 



3 r 



i=l 



1 rll W 



AT|| \ - . 2 



2 

bJl^, (D25) 



- (nx) 



. 1 + A 

Jtt> x 

2 



(D17) X E ^3(™f") 2 -^ 2 +XaH 2 



BjlLf*', (D26) 



■ 9 f 

, rll <! , wy 

Ml = -ig 3 V> - 2wq?7 + — p-^— + 

m 1 — q 



+ iuY,{A l -2 Vl m^)L { r\ (D27) 



Mr 



E 



ATI , 1 + A 

-v 3 m^ H — x 



K 3 (m^) 2 - K 2 q 2 + Xa H 2 ] B±} L™, (D19) 



M 



II _ „ I JV _ r\ 



g\p N -p'\ -(<J a + W )q z -K a q /l + 
3 



1 7 J 

uj p m 
q (m 1 — q) 



^U-2,imf"Uf \ (D28) 
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Ml = ujqWo 



+ iLoY,(A l -2v 1 m^)L[ nx \ (D29) 



Mg 3 = 



N\\ 



Wo - iwif ) flf + 



+ - [vKo + K q 2 - iprnf 11 + iw 7 f 



} 



(nx) 



(D32) 



A/. 



i=l 



+ - \W + K q 2 - ^>mf 11 + iurif 



}L[ VX \ (D30) 



+1 \w a + Ko q 2 ~ 4>mp + iuryf] } L\ vz \ (D31) 



Note that gravity g has been incorporated into the dis- 
persion relation by adding the hydrostatic pressure term 



3 P 



N 



p 1 ] to (Eq. (|D28|) ). 



By setting ^i = 1*2 = ^3, and setting to zero quanti- 
ties Ki, K2, K3, A, 71, 7^,, and Xa, specific to nematic, 
and neglecting curvature contributions by setting to zero 
kq and ipj the dispersion relation is reduced to the well 
studied form for the case of isotropic liquids 0, H, H, i| . 
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